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We study microscopic and macroscopic dynamical proper- 
ties of a frustrated lattice gas showing the violation of Stokes- 
Einstein law. The glassy behaviors are analyzed and related 
with experimental results in glass former systems. 

In simple liquids the connections between the times 
of macroscopic relaxations and the properties of micro- 
scopic particles motion are elucidated by the Stokes- 
Einstein relation which links shear viscosity and parti- 
cle diffusivity. Recent experimental evidence shows how 
such a relation is violated in supercooled glass forming 
liquids. In these systems complex dynamical behaviors, 
which are distributed on many time scales, are supposed 
to be linked in a non trivial way to atomic motions of the 
liquid components [fOP]. 

We face these questions in the context of a microscopic 
model recently introduced to describe some general fea- 
tures of glasses [||. This model, which bridges Spin 
Glasses (SG) |jj and site Frustrated Percolation ||, con- 
sists in a lattice gas model in which each particle contains 
an internal degree of freedom characterized by a spin vari- 
ables. The spins interact via quenched ferromagnetic and 
antiferromagnetic interactions randomly distributed like 
in the SG model. The model is described by the following 
dilute SG Hamiltonian: 

(3H = - J )^(eijSjSj - l)n l n J - p^ m (1) 

(«) i 

where rij = 0,1 are occupancy variables which have an 
internal Ising degree of freedom Si H, the = ±1 are 
quenched random lattice interactions and /i is an adi- 
mensional chemical potential, which plays the role of the 
inverse of the temperature 1/T. 

When all the sites are occupied (rij = 1), i.e. if ^ ^ oo, 
this model is the standard ± J Ising Spin Glass. In the 
limit, J — > oo, it describes a frustrated lattice gas where 
two nearest neighbor particles can be occupied only if 
their internal variables satisfy the constraint eijSiSj = 1. 
Since in a frustrated interaction loop the spins cannot 
satisfy all the couplings, in this model particle configura- 
tions in which a frustrated loop is fully occupied are 
not allowed . 

In a liquid the internal degrees of freedom Si may be 
associated for example to the internal orientation of par- 
ticles, but, more generally, they represent microscopic 
quantities which probe local effects of "frustration" . Ge- 
ometrical hindrance in the liquid implies the existence of 
loops which cannot be fully occupied, in correspondence 
with the frustrated loops in the model (see ||). 



In mean field |7|] the static properties of the model 
show two transitions: one at a value /i c which signals the 
appearance of metastable states and a second thermody- 
namic transition at a higher value /iq (corresponding to a 
lower temperature). These properties recall those of the 
p-Spin Glass models ||]. These models have received re- 
newed attention, as their dynamical behaviour in mean 
field is described by the same equation as those of the 
mode coupling theory for simple liquids Q. Moreover 
the dynamical transition coincide with the first transition 
from the stable to the metastable state. The similarity in 
the static properties of the two models suggest that also 
the dynamics may be the same. Therefore we expect for 
the lattice gas model (|l|) in mean field a dynamical tran- 
sition at /i c with a diffusion coefficient vanishing with a 
power law. 

We have studied the model ([!]) in d = 3 dimensions in 
the limit J — > oo by means of Monte Carlo methods. We 
use a standard Monte Carlo Dynamics, in which particles 
diffuse and spins are updated with Metropolis spin flip. 
Some results on the static and dynamic properties of the 
model have been previously reported ||. Here we focus 
mainly on the relation between the diffusion coefficient 
and the characteristic relaxation time. 

Our data are consistent with the following picture: we 
find a region /i < (x p (high temperature) where the model 
behaves as a normal fluid characterized by single expo- 
nential relaxation and normal diffusion. The value of \i v 
is numerically consistent with the percolation threshold 
for the particles system. For values /i p < [i < the 
density density time dependent autocorrelation function 
exhibits a two step relaxation, corresponding to the a 
and (3 relaxation times. The value fio corresponds to the 
spin glass transition characterized by the divergence of 
the non linear susceptibility and where moreover the dif- 
fusion coefficient, D(n), vanishes In this region D{y) 
exhibits a net cusp at jj,* (/i p < /i* < fj,o). Both the 
diffusion coefficient and the a relaxation time exhibit a 
power law behavior for (i < /x* with a sharp crossover 
to an Arrhenius or Vogel-Fulcher law for /i > //*. This 
crossover can be interpreted in the following way ||. In 
mean field the dynamical transition occurs at the same 
temperature where the static exhibits a transition from 
a stable to a metastable state. Once in the metastable 
state the system is trapped forever due to the infinite 
lifetime of the mean field metastable state. In finite di- 
mension the metastable state has a finite lifetime, and 
the system does not reach a full arrest. 

The value of fj, p is located roughly where the j3 relax- 
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ation and the a relaxation time starts to separate and 
where the relation between the a relaxation time and 
the diffusion coefficient violates the Stokes- Einstein law, 
as observed in real structural glasses |lC|-]l4|] . 

We have considered a cubic lattice of linear sizes 
L = 8, 12, 16 with periodical boundary conditions, fix- 
ing J = 10 and using p as a variable. Our results do not 
change for J = 10000, showing that we are effectively 
close to the value of J = oc. As stated, the dynamic 
in our model is a standard Monte Carlo (MC) dynam- 
ics, in which particle diffuse and spins are updated with 
Metropolis spin flip. The system is equilibrated after suc- 
cessive thermalization at higher and higher values of p for 
about 10 6 MC steps at each fixed value of external pa- 
rameters. The measures are then taken up to about 10 8 
MC steps, for a given random configuration of the cou- 
plings €ij. We are aware that in the deep glassy region 
(above p, ~ 5) our numerical results may be only indica- 
tive due to the required very long simulation times, but 
this do not change the overall picture described above. 

We have calculated the density-density time dependent 
autocorrelation functions: 



on the details of the fit, such as for example the region 



c k {t) = ( Pk (t)p- k m/( Pk p- k ) 



(2) 



Here p k is the Fourier transform on the lattice of the 
density p k (t) = ^ X)r=o ^(M) cos(q fe r), where p{r,t) is 
the particle density, at time t and at a distance r from a 
median plane in our cubic lattice of size L, and q k = ^k, 
with k G {1, 2, L/2}. Times are measured in such a 
way that t = 1 corresponds on average to the time to 
update once all the degree of freedom in the system. Wc 
have also calculate the the square magnetization time de- 
pendent autocorrelation function and found that behaves 
in a similar way. 

As anticipated above, the system at low p (i.e. low 
density) behaves as normal liquid, with exponentially de- 
caying time correlation functions (see Fig. Ill): C k (t) ~ 
cxp(— t/ro(k)). The characteristic time To(k), exclud- 
ing finite size effects, is inversely proportional to k 2 , 
T (k) ~ l/D Q k 2 + Too . 

For larger values of p (larger density) even if the short 
time decay (t < 1) is still exponential, the long time 
relaxation (t >> 1) may be reasonably fitted by the 
Kohlrausch- Williams- Watt stretched exponential form 
(see Fig. |): 



C fe (i) ~ BeM-(t/r(k)f) 



(3) 



with (3 function of p (see Fig. Q), as typically occurs in 
SG |16|. This implies that, the MCT prediction of time- 
temperature superposition relation jl| is verified only in 
very narrow density intervals. 

Close to p p ~ 0.75 (p p ~ 0.38), we observe the separa- 
tion of short times (/3-process) and long times (a-process) 
relaxation. For an estimate of the long time relaxation, 
the a-relaxation, in principle we could use r(k) from the 
fit of cq. (3). However, due to the sensitivity of r(fc) 



chosen for the fit, following Ref. 
laxation T a (k; p) as the time sue 



|] we define the a re- 

h that C k {r a ) = 10" 1 . 



For p < p* ~ 2, r Q (fc] 
behavior (see Fig. ||) JlJ 



is well fitted with a power law 



T a (k; p) = A k {p c - p) 7 + B k 



(4) 



with 7 = 6.8 and p c = 5.6 almost independent of k, 
(A k = 7000 and B k = 0.5 for k = L/4). Also the time 
To (A:), of the initial exponential like decay, follows the 
same power law dependence up to p p , and then saturates 
to a finite value exponentially with p. 

For p > p* , while To(k; p) saturates, r a (k; p) is rea- 
sonable fitted by an Arrhenius behavior divergence (see 
Fig. |): 



Y {k\p) = a k exp(/x//i c 



(5) 



with a k = 0.7 and p a — 0.4 for k = L/4. 

The crossover from a power law behavior forecasted by 
MCT or Spin Glass mean field theory to an Arrhenius (or 
Vogel-Tamman-Fulcher in fragile liquids) behavior for a 
relaxation times, is typically observed in real structural 
glasses too (2|. 

It is interesting to compare the anomalies of the au- 
tocorrelation functions described above, to microscopic 
particle diffusion. Here frustration may have a strong ef- 
fect on particles. We have calculated the particle mean 
square displacement §, R 2 {t) = Ei(n(*) - ^(0)) 2 ). 
For p < p p we find a Brownian typical linear time be- 
havior. For p > ji p , in the intermediate time region, 
we observe a subdiffusive region with an inflection which 
becomes more evident as p increases (see Fig.^J). 

The linear asymptotic behavior of R 2 (t), defines the 
diffusivity D. This shows an apparent cusp at p* ~ 2 
(p* ~ 0.5) and an abrupt change in behavior, as shown 
in Fig.|| Below p* , it is possible to fit D{p) with the 
power law given in eq.Q) with the same 7 and p c (see 
Fig. |^) found studying density relaxation times. 

The value p c from the power law fit corresponds to the 
characteristic temperature T c of Mode Coupling Theory 
jl| , or to the "dynamic transition" of the mean-field the- 
ory of p-Spin Glasses Q,^]. Above p* the best fit for D 
is obtained using a Voghel-Tamman-Fulcher (see Fig. ||): 



D 1 (p) = A D exp(d/(p- 
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)) 



(G) 



with An = 17, p\ = 11.4 and d = 0.3, for the system 
of size 16 3 . However an Arrhenius fit (i.e. p\ 00), as 
that found for r a of a relaxation, is just slightly worse: 



D Hm) = A D exp(p/p D ) 



(7) 



where A' D = 2.72 and pd = 1-1. A crossover from power 
law to Arrhenius (or Voghel-Tamman-Fulcher) behavior 
is also observed in real experiments [||. 

In our model the density p* corresponding to p* , where 
the Arrhenius region sets in, is signaled by a local mini- 
mum in the static structure factor S(k; p), as a function of 
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p. Our observations indicate that the diffusivity D goes 
to zero in the region where the static SG transition, po, 
should be located as signaled by the divergent non linear 
susceptibility xsg- However it is difficult to numerically 
establish where exactly po is located (po > 5.5) (|]. As 
in real glasses, however there is no divergence of density 
fluctuations. 

We compare now the diffusion coefficient D and the 
a relaxation time r a . In a normal liquid the diffusion 
coefficient and the viscosity, 77, are related by the Stokes- 
Einstein relation (SE): r\D = CT where T is the temper- 
ature in Kelvin and C a constant [ fl9f . In our system it is 
not possible to directly define a viscosity, but due to the 
Maxwell relation 77 is proportional to the time scale of 
the asymptotic relaxation r a . We find that the product 
T a D is not a constant varying the potential p, as shown 
in in Fig. |^. Studying T a (p)D(p), three regions emerge, 
separated by the values of p p and p*. A similar separa- 
tion in three different regions is found in real experiments 
on colloidal suspensions 0, where, plotting the product 
r]D as a function of the volume fraction, a behaviour is 
found close to that depicted in the inset of our Fig. [l| 

It has been suggested that the departure from the 
Stokes Einstein relation, in glass forming liquids p0[ , 
could be described by the following relation: D^ 1 = Krf. 
To check if such relation holds in our model we have plot 
(see Fig. f|) 1/D as a function of r Q (fc = L/4). By as- 
suming 

D- 1 = Kri (8) 

we find for p < p p £ = 1 (as in usual SE) , and for p > p* 
£ ~ 0.3. This last value of £ for p > p* is consistent with 
eq. (H) and eq. @) and £ = p a /pD 

Experiments in glass forming liquids jl0]-[l3| indicates 
that the exponent is £ < I, ranging in a broad interval 
depending on the system (for instance in o-terphenyl us- 
ing the rotational diffusion coefficient is found £ = 0.28 
0, in PMMA f = 0.69 @). It is striking that the data 
on o-terphenyl [ |10| , not only exhibit the same exponent 
£ found here, but the rotational diffusion coefficient also 
shows a cusp as function of I/T and the appearance of 
several regions, analogous to those depicted respectively 
in Fig.2 and Fig. §. 

In conclusion we have studied the macroscopic relax- 
ation and microscopic diffusive properties of the frus- 
trated Ising lattice gas introduced in ||. Many connec- 
tions have appeared with the physics of real glass form- 
ing liquids ranging from the anomalies in the density re- 
laxations, to those of diffusivity, to violation of Stokes- 
Einstein law [0J|]. 

The authors are grateful to Dino Leporini for useful 
discussions. This research has been partially supported 
by a grant of CNR. 
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FIG. 1. Fourier transformed density correlation function, 
Cfc(t), for k = L/4, as a function of time t for several values 
of the chemical potential p (p = —2.0,3.0,5.0), in a system 
of linear size L = 16 and coupling constant J = 10. Su- 
perimposed are the short time exponential and the long time 
stretched exponential fits (see text). Inset: The exponent /3 of 
the stretched exponentials of the long time fits of density cor- 
relation functions reported in the main frame, as a function 

of /!. 

FIG. 2. The short time exponential /3-relaxation time 
ro(k) (circle), and the a-relaxation = r a /ln(10) (square), 
of density correlation function (for k = L/4), as a function of 
p, for a system of linear size L = 16 and J = 10. Continuous 
curves are power law and Arrhenius fit described in the text. 
The point where the low p behaviors fails is located around 
p v ~ 0.75. Inset: The inverse of diffusivity, D(p), as a func- 
tion of p for the same system. Superimposed are the power 
law and Voghel-Tamman-Fulcher fits described in the text. 
The cusp in D(p) individuates a characteristic value p* ~ 2.0 
(p* < po, where po is the SG transition). 

FIG. 3. Particles mean square displacement, R(t) 2 , as a 
function of time in a system of linear size L = 16, for J = 10, 
and for p = 0.088, 0.271, 0.440, 0.581, 0.674 (higher curves cor- 
respond to lower densities) . 

FIG. 4. The diffusivity, D, as a function of the a-relaxation 
time, r a . As in real experiments three different regions ap- 
pear. The superimposed curves are fit with the generalized 
fractional Stokes-Einstein law D^ 1 ~ tJ. Inset: The prod- 
uct of diffusivity, D(p), and of a-relaxation time, r a (p), as a 
function of density p. 
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We study microscopic and macroscopic dynamical proper- 
ties of a frustrated lattice gas showing the violation of Stokes- 
Einstein law. The glassy behaviors are analyzed and related 
with experimental results in glass former systems. 

In simple liquids the connections between the times 
of macroscopic relaxations and the properties of micro- 
scopic particles motion are elucidated by the Stokes- 
Einstein relation which links shear viscosity and parti- 
cle diffusivity. Recent experimental evidence shows how 
such a relation is violated in supercooled glass forming 
liquids. In these systems complex dynamical behaviors, 
which are distributed on many time scales, are supposed 
to be linked in a non trivial way to atomic motions of the 
liquid components [1,2]. 

We face these questions in the context of a microscopic 
model recently introduced to describe some general fea- 
tures of glasses [3]. This model, which bridges Spin 
Glasses (SG) [4] and site Frustrated Percolation [5], con- 
sists in a lattice gas model in which each particle contains 
an internal degree of freedom characterized by a spin vari- 
ables. The spins interact via quenched ferromagnetic and 
antiferromagnetic interactions randomly distributed like 
in the SG model. The model is described by the following 
dilute SG Hamiltonian: 

(3% = -j'^2(e ij S i Sj - YjniUj - H^rii (1) 

where n, = 0,1 are occupancy variables which have an 
internal Ising degree of freedom 5, ± 1, the e,j = ±1 are 
quenched random lattice interactions and /j is an adi- 
mensional chemical potential, which plays the role of the 
inverse of the temperature 1/T. 

When all the sites are occupied (n, = 1), i.e. if /j — > oo, 
this model is the standard ±J Ising Spin Glass. In the 
limit, J — > oo, it describes a frustrated lattice gas where 
two nearest neighbor particles can be occupied only if 
their internal variables satisfy the constraint e,j SiSj = 1. 
Since in a frustrated interaction loop the spins cannot 
satisfy all the couplings, in this model particle configura- 
tions in which a frustrated loop [4] is fully occupied are 
not allowed . 

In a liquid the internal degrees of freedom 5, may be 
associated for example to the internal orientation of par- 
ticles, but, more generally, they represent microscopic 
quantities which probe local effects of "frustration" . Ge- 
ometrical hindrance in the liquid implies the existence of 
loops which cannot be fully occupied, in correspondence 
with the frustrated loops in the model (see [6]). 



In mean field [7] the static properties of the model 
show two transitions: one at a value /i c which signals the 
appearance of metastable states and a second thermody- 
namic transition at a higher value liq (corresponding to a 
lower temperature). These properties recall those of the 
p-Spin Glass models [8]. These models have received re- 
newed attention, as their dynamical behaviour in mean 
field is described by the same equation as those of the 
mode coupling theory for simple liquids [1]. Moreover 
the dynamical transition coincide with the first transition 
from the stable to the metastable state. The similarity in 
the static properties of the two models suggest that also 
the dynamics may be the same. Therefore we expect for 
the lattice gas model (1) in mean field a dynamical tran- 
sition at /i c with a diffusion coefficient vanishing with a 
power law. 

We have studied the model (1) in d = 3 dimensions in 
the limit J — > oo by means of Monte Carlo methods. We 
use a standard Monte Carlo Dynamics, in which particles 
diffuse and spins are updated with Metropolis spin flip. 
Some results on the static and dynamic properties of the 
model have been previously reported [3]. Here we focus 
mainly on the relation between the diffusion coefficient 
and the characteristic relaxation time. 

Our data are consistent with the following picture: we 
find a region fx < /j p (high temperature) where the model 
behaves as a normal fluid characterized by single expo- 
nential relaxation and normal diffusion. The value of fi p 
is numerically consistent with the percolation threshold 
for the particles system. For values /j p < /j < /Jo the 
density density time dependent autocorrelation function 
exhibits a two step relaxation, corresponding to the a 
and fi relaxation times. The value /io corresponds to the 
spin glass transition characterized by the divergence of 
the non linear susceptibility and where moreover the dif- 
fusion coefficient, D(p), vanishes [3]. In this region D([i) 
exhibits a net cusp at /j* (/j p < /j* < /io). Both the 
diffusion coefficient and the a relaxation time exhibit a 
power law behavior for /i < /t* with a sharp crossover 
to an Arrhenius or Vogel-Fulcher law for /i > /i*. This 
crossover can be interpreted in the following way [9]. In 
mean field the dynamical transition occurs at the same 
temperature where the static exhibits a transition from 
a stable to a metastable state. Once in the metastable 
state the system is trapped forever due to the infinite 
lifetime of the mean field metastable state. In finite di- 
mension the metastable state has a finite lifetime, and 
the system does not reach a full arrest. 

The value of /j p is located roughly where the f3 relax- 
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ation and the a relaxation time starts to separate and 
where the relation between the a relaxation time and 
the diffusion coefficient violates the Stokes- Einstein law, 
as observed in real structural glasses [10-14]. 

We have considered a cubic lattice of linear sizes 
L = 8, 12, 16 with periodical boundary conditions, fix- 
ing J = 10 and using p as a variable. Our results do not 
change for J = 10000, showing that we are effectively 
close to the value of J = oo. As stated, the dynamic 
in our model is a standard Monte Carlo (MC) dynam- 
ics, in which particle diffuse and spins are updated with 
Metropolis spin flip. The system is equilibrated after suc- 
cessive thermalization at higher and higher values of p for 
about 10 6 MC steps at each fixed value of external pa- 
rameters. The measures are then taken up to about 10 8 
MC steps, for a given random configuration of the cou- 
plings €ij. We are aware that in the deep glassy region 
(above p ~ 5) our numerical results may be only indica- 
tive due to the required very long simulation times, but 
this do not change the overall picture described above. 

We have calculated the density-density time dependent 
autocorrelation functions: 

C k {t) = {p k (t)p-k(0))/{pkP-k) (2) 

Here pk is the Fourier transform on the lattice of the 
density p k (t) = ± Ylr=o P( r > *) cos(q k r), where p(r,t) is 
the particle density, at time t and at a distance r from a 
median plane in our cubic lattice of size L, and q k = ^k, 
with k £ {1, 2, L/2}. Times are measured in such a 
way that t = 1 corresponds on average to the time to 
update once all the degree of freedom in the system. We 
have also calculate the the square magnetization time de- 
pendent autocorrelation function and found that behaves 
in a similar way. 

As anticipated above, the system at low p (i.e. low 
density) behaves as normal liquid, with exponentially de- 
caying time correlation functions (see Fig. 1): Ck(t) ~ 
exp(— i/ro(fc)). The characteristic time ro(fc), exclud- 
ing finite size effects, is inversely proportional to k 2 , 
T (k) ~ 1/D k 2 + 

For larger values of p (larger density) even if the short 
time decay (t < 1) is still exponential, the long time 
relaxation (t >> 1) may be reasonably fitted by the 
Kohlrausch- Williams- Watt stretched exponential form 
(see Fig. 1): 

C k (t) ~Bexp(-(t/T(k)f) (3) 

with (3 function of p (see Fig. 1), as typically occurs in 
SG [16]. This implies that, the MCT prediction of time- 
temperature superposition relation [1] is verified only in 
very narrow density intervals. 

Close to p p ~ 0.75 (p p ~ 0.38), we observe the separa- 
tion of short times (/3-process) and long times (a-process) 
relaxation. For an estimate of the long time relaxation, 
the a-relaxation, in principle we could use r(k) from the 
fit of eq. (3). However, due to the sensitivity of r{k) 



on the details of the fit, such as for example the region 
chosen for the fit, following Ref. [15] we define the a re- 
laxation T a (k;p) as the time such that C k (T a ) = 10 -1 . 
For p < p* ~ 2, T a (k) is well fitted with a power law 
behavior (see Fig. 2) [17]: 

T a (k;p) = A k (p c - p)~' r +B k (4) 

with 7 = 6.8 and p c = 5.6 almost independent of k, 
(A k = 7000 and B k = 0.5 for k = L/4). Also the time 
To(fc), of the initial exponential like decay, follows the 
same power law dependence up to p p , and then saturates 
to a finite value exponentially with p. 

For p > p*, while To(k; p) saturates, r a (k;p) is rea- 
sonable fitted by an Arrhenius behavior divergence (see 
Fig. 2): 

T a (k;p)=a k exp(p/p a ) (5) 

with a k = 0.7 and p a = 0.4 for k = L/4. 

The crossover from a power law behavior forecasted by 
MCT or Spin Glass mean field theory to an Arrhenius (or 
Vogel-Tamman-Fulcher in fragile liquids) behavior for a 
relaxation times, is typically observed in real structural 
glasses too [2]. 

It is interesting to compare the anomalies of the au- 
tocorrelation functions described above, to microscopic 
particle diffusion. Here frustration may have a strong ef- 
fect on particles. We have calculated the particle mean 
square displacement [3], R 2 (t) = (± - r;(0)) 2 ). 

For p < p p we find a Brownian typical linear time be- 
havior. For p > p p , in the intermediate time region, 
we observe a subdiffusive region with an inflection which 
becomes more evident as p increases (see Fig. 3). 

The linear asymptotic behavior of R 2 (t), defines the 
diffusivity D. This shows an apparent cusp at p* ~ 2 
(p* ~ 0.5) and an abrupt change in behavior, as shown 
in Fig. 2. Below p* , it is possible to fit D(p) with the 
power law given in eq.(4) with the same 7 and p c (see 
Fig. 2) found studying density relaxation times. 

The value p c from the power law fit corresponds to the 
characteristic temperature T c of Mode Coupling Theory 
[1], or to the "dynamic transition" of the mean-field the- 
ory of p-Spin Glasses [4,8]. Above p* the best fit for D 
is obtained using a Voghel-Tamman-Fulcher (see Fig. 2): 

D- 1 (p) = A D eMd/ip- 1 - P^ 1 )) (6) 

with Ad = 17, pi = 11.4 and d = 0.3, for the system 
of size 16 3 . However an Arrhenius fit (i.e. p\ — > 00), as 
that found for r a of a relaxation, is just slightly worse: 

D-^p) =A' D exp(n/n D ) (7) 

where A' D =2.72 and pu = 1-1. A crossover from power 
law to Arrhenius (or Voghel-Tamman-Fulcher) behavior 
is also observed in real experiments [2]. 

In our model the density p* corresponding to p* , where 
the Arrhenius region sets in, is signaled by a local mini- 
mum in the static structure factor S(k; p), as a function of 



2 



p. Our observations indicate that the diffusivity D goes 
to zero in the region where the static SG transition, /to, 
should be located as signaled by the divergent non linear 
susceptibility xsg- However it is difficult to numerically 
establish where exactly /Jo is located (/Jo > 5.5) [3]. As 
in real glasses, however there is no divergence of density 
fluctuations. 

We compare now the diffusion coefficient D and the 
a relaxation time r a . In a normal liquid the diffusion 
coefficient and the viscosity, rj, are related by the Stokes- 
Einstein relation (SE): rjD = CT where T is the temper- 
ature in Kelvin and C a constant [19]. In our system it is 
not possible to directly define a viscosity, but due to the 
Maxwell relation rj is proportional to the time scale of 
the asymptotic relaxation r a . We find that the product 
r a D is not a constant varying the potential /j, as shown 
in in Fig. 4. Studying T a (n)D(n), three regions emerge, 
separated by the values of [i p and fi* . A similar separa- 
tion in three different regions is found in real experiments 
on colloidal suspensions [14], where, plotting the product 
r/D as a function of the volume fraction, a behaviour is 
found close to that depicted in the inset of our Fig. 4. 

It has been suggested that the departure from the 
Stokes Einstein relation, in glass forming liquids [10], 
could be described by the following relation: D~ l =Krf-. 
To check if such relation holds in our model we have plot 
(see Fig. 4) 1/D as a function of T a (k = L/A). By as- 
suming 

D- 1 = Kri (8) 

we find for /j < /j p £ = 1 (as in usual SE) , and for /j > /j* 
£ ~ 0.3. This last value of £ for /j > /j* is consistent with 
eq. (5) and eq. (7) and £ = /i a //i_D 

Experiments in glass forming liquids [10-13] indicates 
that the exponent is £ < 1, ranging in a broad interval 
depending on the system (for instance in o-terphenyl us- 
ing the rotational diffusion coefficient is found £ = 0.28 
[10], in PMMA £ = 0.69 [13]). It is striking that the data 
on o-terphenyl [10], not only exhibit the same exponent 
£ found here, but the rotational diffusion coefficient also 
shows a cusp as function of 1/T and the appearance of 
several regions, analogous to those depicted respectively 
in Fig. 2 and Fig. 4. 

In conclusion we have studied the macroscopic relax- 
ation and microscopic diffusive properties of the frus- 
trated Ising lattice gas introduced in [3]. Many connec- 
tions have appeared with the physics of real glass form- 
ing liquids ranging from the anomalies in the density re- 
laxations, to those of diffusivity, to violation of Stokes- 
Einstein law [1,2]. 

The authors are grateful to Dino Leporini for useful 
discussions. This research has been partially supported 
by a grant of CNR. 
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FIG. 1. Fourier transformed density correlation function, 
Ck(t), for k = L/4, as a function of time t for several values 
of the chemical potential 11(11 = —2.0,3.0,5.0), in a system 
of linear size L = 16 and coupling constant J = 10. Su- 
perimposed are the short time exponential and the long time 
stretched exponential fits (see text). Inset: The exponent /3 of 
the stretched exponentials of the long time fits of density cor- 
relation functions reported in the main frame, as a function 
of jl. 



FIG. 3. Particles mean square displacement, R(t) 2 , as a 
function of time in a system of linear size L = 16, for J = 10, 
and for p = 0.088, 0.271, 0.440, 0.581, 0.674 (higher curves cor- 
respond to lower densities). 
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FIG. 2. The short time exponential /3-relaxation time To{k) 
(circle), and the a-relaxation = r a / ln(10) (square), of den- 
sity correlation function (for k = L/4), as a function of p, for 
a system of linear size L = 16 and J = 10. Continuous 
curves are power law and Arrhenius fit described in the text. 
The point where the low fj, behaviors fails is located around 
fi p ~ 0.75. Inset: The inverse of diffusivity, D(ji), as a func- 
tion of fi for the same system. Superimposed are the power 
law and Voghel-Tamman-Fulcher fits described in the text. 
The cusp in D(ji) individuates a characteristic value ~ 2.0 
(p* < fio, where /io is the SG transition). 
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FIG. 4. The diffusivity, D, as a function of the a-relaxation 
time, To,. As in real experiments three different regions ap- 
pear. The superimposed curves are fit with the generalized 
fractional Stokes-Einstein law Z) _1 ~ r|. Inset: The prod- 
uct of diffusivity, D(p), and of a-relaxation time, r a (p), as a 
function of density p. 
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